Cost-effective Whole Exome Sequencing discovers pathogenic variant causing Neurofibromatosis type 1 in a family from Jammu and Kashmir, India.

Neurofibromatosis type 1 (NF1) is a multisystemic hereditary disorder associated with an increased risk of benign and malignant tumor formation predominantly on the skin, bone, and peripheral nervous system. It has been reported that out of all the NF1 cases, more than 95% cases develop the disease due to heterozygous loss-of-function variants in Neurofibromin (NF1) gene. However, identification of NF1 causative variants by presently recommended method of gene-targeted Sanger sequencing is challenging and cost-intensive due to the large size of the NF1gene with 60 exons spanning about 350 kb. Further, conducting the genetic studies is difficult in low resource regions and among families with the limited financial capabilities, restricting them from availing diagnostic as well as proper disease management measures. Here, we studied a three-generation family from Jammu and Kashmir state in India, with multiple affected family members showing clinical indications of NF1. We combinedly used two applications, Whole Exome Sequencing (WES) and Sanger sequencing, for this study and discovered a nonsense variant NM_000267.3:c.2041C>T (NP_000258.1:p.Arg681Ter*) in exon 18 of NF1 gene in a cost effective manner. In silico analyses further substantiated the pathogenicity of this novel variant. The study also emphasized on the role of Next Generation Sequencing (NGS) as a cost-effective method for the discovery of pathogenic variants in disorders with known phenotypes found in large sized candidate genes. The current study is the first study based on the genetic characterization of NF1 from Jammu and Kashmir-India, highlighting the importance of the described methodology adopted for the identification and understanding of the disease in low resource region. The early diagnosis of genetic disorders would open the door to appropriate genetic counseling, reducing the disease burden in the affected families and the general population at large.

Genetic interaction between F-box motif encoding YDR131C and retrograde signaling-related RTG1 regulates the stress response and apoptosis in Saccharomyces cerevisiae.

The retrograde signaling pathway is well conserved from yeast to humans, which regulates cell adaptation during stress conditions and prevents cell death. One of its components, RTG1 encoded Rtg1p in association with Rtg3p communicates between mitochondria, nucleus, and peroxisome during stress for adaptation, by regulation of transcription. The F-box motif protein encoded by YDR131C  constitutes a part of SCF Ydr131c -E3 ligase complex, with unknown function; however, it is known that retrograde signaling is modulated by the E3 ligase complex. This study reports epistasis interaction between YDR131C and RTG1, which regulates cell growth, response to genotoxic stress, decreased apoptosis, resistance to petite mutation, and cell wall integrity. The cells of ydr131cΔrtg1Δ genetic background exhibits growth rate improvement however, sensitivity to hydroxyurea, itraconazole antifungal agent and synthetic indoloquinazoline-based alkaloid (8-fluorotryptanthrin, RK64), which disrupts the cell wall integrity in Saccharomyces cerevisiae. The epistatic interaction between YDR131C and RTG1 indicates a link between protein degradation and retrograde signaling pathways.

Review: Understanding Rare Genetic Diseases in Low Resource Regions Like Jammu and Kashmir - India.

Rare diseases (RDs) are the clinical conditions affecting a few percentage of individuals in a general population compared to other diseases. Limited clinical information and a lack of reliable epidemiological data make their timely diagnosis and therapeutic management difficult. Emerging Next-Generation DNA Sequencing technologies have enhanced our horizons on patho-physiological understanding of many of the RDs and ushered us into an era of diagnostic and therapeutic research related to this ignored health challenge. Unfortunately, relevant research is meager in developing countries which lack a reliable estimate of the exact burden of most of the RDs. India is to be considered as the "Pandora's Box of genetic disorders." Owing to its huge population heterogeneity and high inbreeding or endogamy rates, a higher burden of rare recessive genetic diseases is expected and supported by the literature findings that endogamy is highly detrimental to health as it enhances the degree of homozygosity of recessive alleles in the general population. The population of a low resource region Jammu and Kashmir (J&K) - India, is highly inbred. Some of its population groups variably practice consanguinity. In context with the region's typical geographical topography, highly inbred population structure and unique but heterogeneous gene pool, a huge burden of known and uncharacterized genetic disorders is expected. Unfortunately, many suspected cases of genetic disorders remain undiagnosed or misdiagnosed due to lack of appropriate clinical as well as diagnostic resources in the region, causing patients to face a huge psycho-socio-economic crisis and many a time suffer life-long with their ailment. In this review, the major challenges associated with RDs are highlighted in general and an account on the methods that can be adopted for conducting fruitful molecular genetic studies in genetically vulnerable and low resource regions is also provided, with an example of a region like J&K - India.

Next generation sequencing and microbiome's taxonomical characterization of frozen soil of north western Himalayas of Jammu and Kashmir, India

BACKGROUND: Traditionally, microbial genome sequencing has been restrained to the species grown in pure culture. The development of culture-independent techniques over the last decade allows scientists to sequence microbial communities directly from environmental samples. Metagenomics is the study of complex genome by the isolation of DNA of the whole community. Next generation sequencing (NGS) of metagenomic DNA gives information about the microbial and taxonomical characterization of a particular niche. The objective of the present research is to study the microbial and taxonomical characterization of the metagenomic DNA, isolated from the frozen soil sample of a glacier in the north western Himalayas through NGS. RESULTS: The glacier community comprised of 16 phyla with the representation of members belonging to Proteobacteria and Acidobacteria. The number of genes annotated through the Kyoto Encyclopedia of Genes and Genomes (KEGG), GO, Pfam, Clusters of Orthologous Groups of proteins (COGs), and FIG databases were generated by COGNIZER. The annotation of genes assigned in each group from the metagenomics data through COG database and the number of genes annotated in different pathways through KEGG database were reported. CONCLUSION: Results indicate that the glacier soil taken in the present study, harbors taxonomically and metabolically diverse communities. The major bacterial group present in the niche is Proteobacteria followed by Acidobacteria, and Actinobacteria, etc. Different genes were annotated through COG and KEGG databases that integrate genomic, chemical, and systemic functional information.

Genetic variants of DNAH11 and LRFN2 genes and their association with ovarian and breast cancer.

OBJECTIVE: To investigate the association of newly identified genetic variants G>A (rs2285947) of the DNAH11 gene and G>A (rs2494938) of the LRFN2 gene with ovarian and breast cancers in women belonging to Jammu and Kashmir state, where the prevalence of ovarian and breast cancers is remarkably high in the population. METHODS: A candidate gene prospective case-control association study design was adopted, in which 354 cases (219 cases of ovarian cancer and 135 cases of breast cancer) were histopathologically confirmed and 330 healthy controls matched for age and ethnicity were recruited. The details of cases and controls were also recorded in a predesigned pro forma after their written informed consent. Both variants were genotyped by TaqMan allele discrimination assay using real-time polymerase chain reaction. Logistic regression analysis was performed to estimate the corrected odds ratio (OR), confidence interval (CI), and level of significance (P value) for potential confounding factors. RESULTS: The rs2285947 variant of DNAH11 was found to be significantly associated with both ovarian and breast cancers with adjusted ORs of 1.7 (95% CI 1.2-2.4; P=0.004) and 1.70 (95% CI 1.13-2.54; P=0.0009), respectively. However, no significant association of variant rs2494938 of LRFN2 was observed with ovarian cancer (estimated OR 0.9, 95% CI 0.6-1.4; P=0.919) or breast cancer (estimated OR 1.27, 95% CI 0.8-1.9; P=0.216). CONCLUSIONS: The collected data proposed that the variant rs2285947 of DNAH11 gene is a potential risk factor for ovarian and breast cancers in the studied population.

A hydrolase with esterase activity expressed from a fosmid gene bank prepared from DNA of a North West Himalayan glacier frozen soil sample.

Screening of 20,000 clones of a fosmid gene bank, constructed from DNA extracted from North West Himalaya (NWH) glacier soil sample, using functional approach identified 10 esterase/lipase-producing clones. Of these, a clone designated pFG43 with an insert size of 45 kb which produced the highest concentration of enzyme (467.43 U/mg) was sequenced. Clone pFG43 contained 61 open reading frames (ORF) and of these an ORF of 1155 bp designated ME-003, was found to be closely related to a hydrolase from Acidobacteria sps (77% sequence identity and E value = 1e-164) and subsequently identified as a putative cocaine esterase. ORF ME-003 was amplified and sub-cloned using a TA vector system into E. coli (DH5α). The purified recombinant enzyme with a molecular weight of 43 kDa had optimal activity at 40 °C, pH 6 and the highest activity with shorter chain fatty acids than with higher chain length fatty acids. There is insignificant effect of inhibitors on the enzyme activity of ME-003, except PMSF which completely inhibited its activity. ME-003 activity was also inhibited in the presence of copper oxide but remained stable in presence of other metal ions. The enzyme activity was also inhibited in the presence of organic solvents; however, in the presence of 10% isopropanol, 12% of enzymatic activity was retained. Among various detergents, SDS completely inhibited enzymatic activity. The recombinant enzyme also shows enantio-specific activity against the racemic drug intermediates/precursors and exhibited 90% ee against racemic 1-phenyl ethanol and fluoxetine.

Molecular dynamics analysis of the effects of GTP, GDP and benzimidazole derivative on structural dynamics of a cell division protein FtsZ from Mycobacterium tuberculosis.

The prevailing multi-drug resistance in Mycobacterium tuberculosis continues to remain one of the main challenges to combat tuberculosis. Hence, it becomes imperative to focus on novel drug targets. Filamenting temperature-sensitive mutant Z (FtsZ) is an essential cell division protein, a eukaryotic tubulin homologue and a promising drug target. During cytokinesis, FtsZ polymerises in the presence of GTP to form Z-ring and recruits other proteins at this site that eventually lead to the formation of daughter cells. Benzimidazoles were experimentally shown to inhibit Mtb-FtsZ, with one of the benzimidazole derivatives, M1, being reported to have the minimum inhibitory concentration (MIC) value of 3.13 µg/mL. In the present study, mechanism of destabilisation of FtsZ in the presence of M1 was computationally investigated in the presence of its substrate GTP/GDP employing molecular dynamics (MD) simulation analysis, principal component analysis (PCA), molecular mechanics combined with the generalised Born and surface area continuum salvation (MM-GBSA) and density functional theory (DFT). From the analyses, it is proposed that binding of M1 in the inter-domain cleft induces structural changes in the GTP-binding region that affect GTP binding, thus switching the preference of this protein towards depolymerised state and eventually inhibiting the cell division. Hence, this study provides mechanistic insights into the design of novel benzimidazole inhibitors against Mtb-FtsZ. Communicated by Ramaswamy H. Sarma.

In Silico Evaluation of Variable pH on the Binding of Epidermal Growth Factor Receptor Ectodomain to its Ligand Through Molecular Dynamics Simulation in Tumors.

Many aggressive and metastatic cancer cell types show Warburg Effect; therefore, it is a possible adaptation helping cancer cells to rapidly divide and utilize the glycolytic intermediates for biosynthesis of ribose sugars (for nucleotide biosynthesis), fatty acid synthesis (lipids for membrane synthesis), NADPH (cellular currency for reductive biosynthesis) and lactate. This in due course results in decrease of extracellular pH, leading to acidic tumor micro-environment. EGFR is a crucial cell surface signaling receptor implicated in cancer cell survival and progression. This warrants studying the effect of the acidic micro-environmental conditions on the binding of the EGFR cell surface receptor to one of its natural extracellular ligand EGF. We exploited in silico approaches: molecular dynamics simulation at variable pH and MM-GBSA free energy of binding calculation method to evaluate the effect of this change in microenvironmental pH. Through the present study it is reported that at pH 6.6 the EGFR binds to EGF with decreased free energy of binding as compared to pH 7.2.

DNA base excision repair genes variants rs25487 (X-ray repair cross-complementing 1) and rs1052133 (human 8-oxoguanine glycosylase 1) with susceptibility to ovarian cancer in the population of the Jammu region, India.

BACKGROUND: Ovarian cancer is highly prevalent in the population of Jammu, in India; the ovarian cancer ranks third among other types of cancer prevalent in females. However, association studies on ovarian cancer are lacking in this region. We aimed to investigate the disease susceptible variants rs1052133 (human 8-oxoguanine glycosylase 1 [hOGG1]) and rs25487 (X-ray repair cross-complementing 1 [XRCC1]) with ovarian cancer in population of Jammu, India. MATERIALS AND METHODS: The study conducted in the Shri Mata Vaishno Devi University is a 3-year study which included a total of 280 well-characterized samples (130 ovarian cancer cases and 150 healthy controls). hOGG1 and XRCC1 polymorphisms were determined by polymerase chain reaction-based restriction fragment length polymorphism, and these genotyping results were confirmed by Sanger sequencing. Hardy-Weinberg equilibrium for both single-nucleotide polymorphisms (SNPs) was assessed using the Chi-square test. The allele and genotype-specific risks were estimated by odds ratios with 95% confidence intervals. RESULTS: In this preliminary study, SNP rs1052133 showed protection with ovarian cancer (P = 0.042). The SNP rs25487 was not found associated with ovarian cancer (P = 0.271). CONCLUSION: Our results indicate that the G allele of rs1052133 imparts protection to the population whereas variant rs25487 was not associated with ovarian cancer in population from the Jammu region, indicating that larger sample size is needed for further statistical validation. Further, association of other SNPs in these genes should also be carried out as their role cannot be ruled out.

Production and purification of an alkaline lipase from Bacillus sp. for enantioselective resolution of (±)-Ketoprofen butyl ester.

The present study was conducted to purify lipase from indigenous Bacillus subtilis strain Kakrayal_1 (BSK-L) for enantioselective resolution of racemic-ketoprofen. The production of lipase (BSK-L) was optimized using Plackett-Burman and central composite design of response surface methodology (RSM). The optimized media containing olive oil (3.5%), MnSO4 (8 mM), CaCl2 (5 mM), peptone (20 g/l), pH (8), agitation (180 rpm) and temperature (37 °C) resulted in maximum lipase production of 7500 U/g of cell biomass. The lipase was purified using sequential method to an overall purification fold of 13% with 20% recovery, 882 U/mg specific activity and a molecular weight of 45 kDa. Optimal pH and temperature of purified lipase were found to be 8 and 37 °C, respectively. Furthermore, BSK-L displayed good stability with various organic solvents, surfactants and metal ions. K m and V max values of lipase were observed to be 2.2 mM and 6.67  mmoles of product formed/min/mg, respectively. The racemic ketoprofen butyl ester was hydrolyzed using lipase with 49% conversion efficiency and 69% enantiomeric excess (ee) which was superior to the commercially procured lipase (Candida antarctica lipase). Thus, this enzyme could be considered as a promising candidate for the pharmaceutical industry.

In silico evaluation of the resistance of the T790M variant of epidermal growth factor receptor kinase to cancer drug Erlotinib.

Epidermal growth factor receptor kinase is implicated in cancer development due to either overexpression or activation variants in its functional intracellular kinase domain. Threonine to methionine (Thr 790 Met) is one such variant observed commonly in patients showing resistance to kinase inhibitor drug Erlotinib. Two mechanisms for resistance have been proposed (1) steric hindrance and (2) enhanced binding to ATP. In this study, we employed molecular dynamics simulations and studied both the mechanisms. Extensive simulations and free energy of binding analyses has shown that steric hindrance does not explain appropriately the mechanism for resistance against Erlotinib therapy for this variant. It has been observed that conformational switching from an intermediate intrinsically disordered C-helix conformation is required for completion of the kinase's catalytic cycle. Our study substantiates that T790M variant has greater tendency for early transition to this intrinsically disordered C-helix intermediate state. We propose that enhanced catalytic efficiency in addition to enhanced ATP binding explains mechanism of T790M resistance to drug Erlotinib.

Evaluation of alkali and thermotolerant lipase from an indigenous isolated Bacillus strain for detergent formulation

Background: Lipases are used in detergent industries to minimise the use of phosphate-based chemicals in detergent formulations. The use of lipase in household laundry reduces environmental pollution and enhances the ability of detergent to remove tough oil or grease stains. Results: A lipase-producing indigenous Bacillus subtilis strain [accession no. KT985358] was isolated from the foothills of Trikuta mountain in Jammu and Kashmir, India. The lipase (BSK-L) produced by this strain expressed alkali and thermotolerance. Lipase has an optimal activity at pH 8.0 and temperature 37°C, whereas it is stable at pH 6.0­9.0 and showed active lipolytic activity at temperatures 30 to 60°C. Furthermore, lipase activity was found to be stimulated in the presence of the metal ions Mn2+, K+, Zn2+, Fe2+ and Ca2+. This lipase was resistant to surfactants, oxidising agents and commercial detergents, suggesting it as a potential candidate for detergent formulation. BSK-L displayed noticeable capability to remove oil stains when used in different washing solutions containing buffer, lipase and commercial detergent. The maximum olive oil removal percentage obtained was 68% when the optimum detergent concentration (Fena) was 0.3%. The oil removal percentage from olive oil-soiled cotton fabric increased with 40 U/mL of lipase. Conclusions: This BSK-L enzyme has the potential for removing oil stains by developing a pre-soaked solution for detergent formulation and was compatible with surfactants, oxidising agents and commercial detergents.

A variation in PANK2 gene is causing Pantothenate kinase-associated Neurodegeneration in a family from Jammu and Kashmir - India.

Pantothenate kinase-associated neurodegeneration is a rare hereditary neurodegenerative disorder associated with nucleotide variation(s) in mitochondrial human Pantothenate kinase 2 (hPanK2) protein encoding PANK2 gene, and is characterized by symptoms of extra-pyramidal dysfunction and accumulation of non-heme iron predominantly in the basal ganglia of the brain. In this study, we describe a familial case of PKAN from the State of Jammu and Kashmir (J&K), India based on the clinical findings and genetic screening of two affected siblings born to consanguineous normal parents. The patients present with early-onset, progressive extrapyramidal dysfunction, and brain Magnetic Resonance imaging (MRI) suggestive of symmetrical iron deposition in the globus pallidi. Screening the PANK2 gene in the patients as well as their unaffected family members revealed a functional single nucleotide variation, perfectly segregating in the patient's family in an autosomal recessive mode of inheritance. We also provide the results of in-silico analyses, predicting the functional consequence of the identified PANK2 variant.

Biotransformation of methyl-ß-D-glucopyranoside to higher chain alkyl glucosides by cell bound ß-glucosidase of Pichia etchellsii.

In the present study, we have investigated the use of Pichia etchellsii whole cells for synthesis of long-chain alkyl glucosides. Methyl-ß-d-glucopyranoside (MG) was used in reaction with fatty alcohols, n-hexanol, n-octanol, n-decanol and n-dodecanol to synthesize the respective alkyl glucosides. The initial reaction conditions were first optimized at 2.5 ml scale for synthesis of octyl glucoside (OG) and were 8% water content, 100mM MG and 6h of reaction time and this resulted in ≈ 53% yield. A maximum transglucosylation/hydrolysis ratio of 2.79 was obtained at 100mM MG favoring high product yield. Based on the optimized conditions, a reactor was operated at 50 ml level which resulted in ≈ 60% conversion of MG to OG. A simple high performance liquid chromatography method was developed for quantitation of higher chain glucosides using a refractive index detector. A maximum of 27% and 13% yield was obtained for decyl-, and dodecyl-ß-d-glucopyranoside, respectively.

In silico proteomic characterization of human epidermal growth factor receptor 2 (HER-2) for the mapping of high affinity antigenic determinants against breast cancer.

Multiple different approaches are being used to activate the immune system against breast cancer. Vaccine therapy in general follows the principle that injections of various substances ultimately result in the presentation of tumor peptides to the patient's immune system. We proposed a potential in silico DNA vaccine against breast cancer by integrating high affinity T cell (MHC-I and MHC-II) and B cell (continuous and discontinuous) epitopes. The matching of the HLA haplotype and antigen was performed to provide the appropriate peptide epitope suitable for majority of the patients. The immunogenic nature of the antigenic construct was also enhanced by the administration of consensus epitopes. The potency of DNA vaccines depends on the efficient expression and presentation of the encoded antigen of interest and the chances of efficient expression of our antigenic construct in host organism was also verified by in silico approaches. An attempt was made to overcome the limited potency of the DNA vaccine by targeting DNA to professional antigen-presenting cells (APCs). A higher immune response theoretically corresponds to a higher survival rate of patients. Therefore, optimization studies were also employed to enhance the immunogenicity of proposed in silico DNA vaccine.

Comparative analysis of mesothelial invasion by single- and multi-wall carbon nanotubes using computational approach.

This study explored the interactions of carbon nanotubes (CNTs) with mesothelial cells. We carried out molecular dynamic simulation studies on the mesothelial cell invasion by single and multi walled carbon nanotubes. In our simulation experiments SWCNTs showed greater binding affinity with membrane in comparison to MWCNTs with similar aspect ratio.

In silico identification of novel protective VSG antigens expressed by Trypanosoma brucei and an effort for designing a highly immunogenic DNA vaccine using IL-12 as adjuvant.

African trypanosomiasis continues to be a major health problem, with more adults dying from this disease world-wide. As the sequence diversity of Trypanosoma brucei is extreme, with VSGs having 15-25% identity with most other VSGs, hence it displays a huge diversity of adaptations and host specificities. Therefore the need for an improved vaccine has become an international priority. The highly conserved and specific epitopes acting as both CD8+ and CD4+ T-cell epitopes (FLINKKPAL and FTALCTLAA) were predicted from large bunch of VSGs of T. brucei. Besides, some other potential epitopes with very high affinity for MHC I and II molecules were also determined while taking consideration on the most common HLA in the general population which accounts for major ethnicities. The vaccine candidates were found to be effective even for non-african populations as predicted by population coverage analysis. Hence the migrating travelers acting as a spread means of the infection can probably also be treated successfully after injection of such a multiepitopic vaccine. Exploiting the immunoinformatics approaches, we designed a potential vaccine by using the consensus epitopic sequence of 388 VSG proteins of T. brucei and performed in silico cloning of multiepitopic antigenic DNA sequence in pBI-CMV1 vector. Moreover, various techniques like codon adaptation, CpG optimization, removal of self recognized epitopes, use of adjuvant and co-injection with plasmids expressing immune-stimulatory molecules were implemented to enhance the immunogenicity of the proposed in silico vaccine.

In silico designing and optimization of anti-breast cancer antibody mimetic oligopeptide targeting HER-2 in women.

Overexpression of HER-2 is of frequent (20-30%) occurrence in breast cancer. Therapeutic targeting of HER-2 with humanized antibody derived oligopeptide may be a promising approach to the treatment of breast cancer. HER-2 gene is part of a family of genes that play critical roles in regulating transmembrane growth of breast cancer cells. Pertuzumab, a recombinant humanized monoclonal antibody (2C4), binds to extracellular domain II of the HER-2 receptor and inhibits its ability to dimerize with other HER receptors blocking the cell growth, signaling and apoptosis induction. The unique binding pocket on HER-2 for pertuzumab provides an important target domain for creation of new anticancer drugs. In the present work an efficient oligopeptide was designed by our computational method that interacts with pertuzumab binding sites of HER-2. In silico docking study demonstrated the best specific interaction of RASPADREV oligopeptide with the dimerization domain in the HER-2 molecule among various screened oligopeptides. ADMET and SAR properties prove the drug likeness of designed oligopeptide as having value 0.98.

Functional cloning and predictive structural modeling of a novel esterase from Bacillus subtilis strain, RRL 1789.

We have recently reported the purification and characterization of a novel esterase from the Bacillus subtilis strain. In the present study we report the genomic DNA cloning and predictive structural modeling of this novel esterase. Tributyrin- and Rhodamine B-based functional screen of a Bacillus subtilis genomic library led to the identification of a potential lipolytic gene. DNA sequence analysis of the cloned gene showed that it encodes a protein of 489 amino acid residues. Sequence homology search and multiple sequence alignment showed that the protein was highly homologous to known esterases. Secondary structure-driven multiple sequence alignment with the homologous esterase of known three-dimensional structures was performed and a 3D structure model of this enzyme was constructed. Based on the topological organization of the secondary structures, this protein belongs to the alpha/beta hydrolase superfamily. Moreover, the presence of serine in the context of amino acid sequence G/A-X-S-X-G (with X an arbitrary amino acid residue) in the protein indicates that it belong to the class of serine hydrolases of this superfamily.

Application of a multiplex PCR assay for the detection of Shigella, Escherichia coli and Shiga toxin-producing Esch. coli in milk.

A multiplex PCR (mPCR) assay using previously known genetic markers of Shigella, Escherichia coli and Shiga-toxic Esch. coli was standardized. uidA gene was targeted for the common detection of Esch. coli and Shigella, whereas ipaH and stx1 genes were used as markers for the detection of Shigella and shiga-toxin producing strains, respectively. The standardized assays detected the target organism specifically and selectively. The mPCR developed by combining all the three reactions generated specific products. The inclusivity and exclusivity tests depicted the precise specificity of the mPCR assay. Results were interpreted on the basis of the pattern of amplicons generated: amplifications of the ipaH and uidA gene fragments indicated the presence of Shigella spp., amplification of uidA alone revealed the presence of Esch. coli and additional presence of verotoxin gene amplicon indicated verotoxinogenic nature of the strain. Specific patterns of bands were obtained when different strains of Esch. coli and Shigella spp. were subjected to this assay. The reactions, individually as well as in the mPCR, could detect approximately 1 cell per 20-microl PCR assay. The protocols were validated by analyzing the coded samples of full fat milk spiked with different pathogens. In naturally contaminated raw milk samples (n=100), Esch. coli were detected in all samples and verotoxinogenic Esch. coli in 15 samples. Shigella, however, was not detected in any of the samples. When DNA purified from the samples found positive for Shiga-toxic Esch. coli was directly used as template for the mPCR, the results showed agreement with the enrichment based detection. The mPCR assay, standardized in this study, may be used for rapid microbiological evaluation of milk samples. Further, the study emphasizes the need for better hygienic conditions in dairies.